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Abstract 

The problem of target localization with ideal binary detectors is considered 
in one dimensional space. The problem is investigated in both a censored 
and non-censored scheme. In the censored setting, the problem is equivalent 
to estimating the center of a uniform distribution by knowing samples of 
data. It does not admit an MVU estimator according to the previous results 
of Lehmann-Sheffe. However, it is proven that if the radius of detection is 
known and sensor deployment region is very large, both censored and non- 
censored cases will have an MVU estimator among the functions that are 
invariant to Euclidean motion. In addition, it is shown that when the radius 
of detection is not known, the censored case still has an MVU estimator 
whereas in the non-censored case, an MVU estimator does not exist, even 
under the assumption that the estimators are invariant to Euclidean motion. 
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1. Introduction 


The problem of localization of a target using a number of distributed 
measurements has been addressed widely in the electrical engineering lit¬ 
erature nna. The problem originally was considered in radar and sonar 
011] but there is renewed interest in it due to advances in wireless sensor 
networks ISH7I and applications in 4G/ LTE networks p-HO]. The localiza¬ 
tion usually is performed through measurements such as Time Difference of 
Arrival (TD0A)[IIHI3|, Direction of Arrival IDOA) [TH - [TB] . Frequency Dif¬ 
ference of Arrival iFDOAl fTTHTU] . and Radio Signal Strength fRSS) [2U] - [51] . 
The hrst three methods do not £t well within energy and complexity con¬ 
straints of wireless sensor networks especially when the target is not coop¬ 
erative [20l [32|. Examples of applications when the target does not assist in 
evaluating its position are identifying the primary user in cognitive radio [31], 
spectrum cartography [3^ , identify unauthorized users of bandwidth and lo¬ 
calizing jammers positions in the battleheld [36]. Although measurements 
of received signal strength can be exact, it is not practical to transmit a 
non-quantized measurement result to the fusion center (FC), where the de¬ 
cision about the location of a target is reached because unlimited bandwidth 
is needed. Thus, in many papers it is assumed that the sensors make a 
quantized [2B] or binary measurement [2SHMI133 EB] of the received signal 
strength. Binary measurements are preferred because the implementation is 
simple and requires minimum bandwidth. In addition, in localization using 
binary sensors, some papers assume that only the detecting sensors (which 
generate a ’’one” output) will communicate data to FC [HIl EH 133 EH] • This 
censoring scheme saves communication overhead and energy because if the 
sensor deployment region is very large it is expected that most sensors do 
not detect the target EDI- 

From a mathematical point of view, the problem can be considered as 
a point-wise estimation when samples of data are available from a popula¬ 
tion function representing the probability of detection at a specihc distance 
from the target. The case when measurements are noise-free is a limiting 
scenario and its performance can be considered as the lower bound for all 
other localizers with uncertain measurements such as those in the presence 
of noise and fading. Thus, analysis of this scenario is important. In this 
case, the problem is equivalent to estimating the centroid of a uniform dis¬ 
tribution from a limited number of samples. However, it appears that this 
formulation of the problem does not have a Cramer Rao Bound (CRB) be- 
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cause the likelihood function is not well behaved [H]. Therefore, a reasonable 
approach to establish the lower bound for the performance in this scenario 
would be to hnd an MVU estimator. This raises the question if a mini¬ 
mum variance unbiased estimator exists at all. One standard approach to 
this problem is to hrst examine the existence of complete sufficient statistics 
|41j . The relation between completeness and existence of a uniform minimum 
variance unbiased estimator has long been of interest to Statisticians. For 
example Lehmann-Sheffe show that if a complete sufficient statistics exists, 
all estimable parametric functions could be uniformly MVU estimated |42] . 
Unfortunately, in this case a complete sufficient statistics does not exist most 
of the time as we will see in Sections 3.1 and|^ Bahadur in |13] provides a 
converse to the Lehmann-Sheffe theorem, and show that if every estimable 
parametric function admits a UMVUE (uniform minimum variance unbiased 
estimator), then a complete sufficient statistic exists. However, |1H 05] pro¬ 
vide examples to show that even when a complete sufficient statistics does 
not exist, there may still be a chance that an MVU estimator exists for some 
parametric function. In Section 3.1, we will mention some known mathe¬ 


matical results that suggests that an MVU estimator does not exist for this 
problem in the censored setting. However, we provide a constraint on the 
function space which seems necessary in location estimation and prove that 
under that constraint, the MVU estimator exists. We believe that this con¬ 
straint is critical for many applied estimators and should be considered in the 
evaluation of MVU estimators in many other applications. We also consider 
different cases of this problem and discuss the existence of MVU estimators 
for each of them. We believe that this work also may have pedagogical value 
in the context of statistical inference. 

In Section]^ we will formalize the problem and mathematical notions. In 
Section we will consider the problem in a censored scenario. In Section 
we will discuss the results for a non-censored case. Concluding remarks and 
discussion is presented in Section 


2. Problem Formulation 

Assume that a target is located at an unknown location zt in space and 
transmits a signal whose power propagates isotropically and is attenuated 
monotonically as a function of distance from the target. To simplify matters, 
we solve the problem in one dimensional space, an approach that has been 
widely studied in the statistics literature. N sensors, are randomly scattered 
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in a deployment region, G = [Ai, A2] of length L = A2 — Ai. They measure 
the received power and compare it to a threshold, r, to make a binary decision 
about the target presence. Let n be the number of detecting sensors (the 
ones that report decision ’’one” to FC) and p be the number of non-detecting 
sensors (the ones whose decisions are ’’zero”). Also let z-^i and represent 
the locations of the ith detecting and jth non-detecting sensors respectively. 
We represent the set of indices of all detecting sensors by iSd = {1,.., n} and 
the set of indices of all non-detecting sensors by iSnd = {1, Therefore, 

Zd = [^Dil* £ ^d] is a vector containing the locations of all detecting sensors 
and Znd = \j ^ i^nd] is a similar vector for non-detecting sensors. 

We assume the sensors make a noise free decision, which can be considered 
as the limiting case when the measured power is averaged over a sufficiently 
long duration of time. Since the received power is a decreasing function 
of distance from the target, this means that all the detecting sensors are 
located within a detection radius, R, from the target i.e. Vz G Sj:,,z-oi G 
[zt — R, zt + R]- R can be considered known or unknown depending on 
whether or not the propagation model and the transmit power are known to 
FC; we will investigate each case separately. We also assume that at least 
one sensor is a detecting sensor {n > 1) and G is sufficiently large such that 
[zt — 3R, zt + 3R] C [^ 1 ,^ 2 ]. Hence, the results of this work are also valid 
when G is unlimited (entire x axis). In addition we assume that at least one 
sensor detects the target. 

We will consider the problem both in a censored and non-censored scheme. 
In a censored scenario presented in Section we assume the sensors are con- 
hgured such that only the detecting sensors report their locations to the FC. 
Whereas, in the non-censored scenario discussed in Section both detecting 
and non-detecting sensors inform FC of their locations and decisions. 
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3. On the Existence of MVU Estimator in the Censored Scheme 

In this scenario, only the locations of detecting sensors are available to 
FC. Therefore, the localization problem is considered as a special instance 
of point estimation and can be viewed as estimating the center of a uniform 
distribution by knowing samples of data from the uniform distribution. 

We divide the problem into the cases in which the radius of detection, 
R, is known to FC and when it is unknown and we investigate each one 
separately. 


3.1. Known Detection Radius 

In this case, R is known and the goal is to estimate zt from n samples 
of a uniform distribution over [zt — R, + R\. It is known in the statistics 
literature that ^tiu) which zdl = niinjg^j^ zj^ and z-qu = maxjg^p z^ 
is a minimal sufficient statistics for estimation of zt gsiiiz]. The intuition 
behind is that if the minimum and maximum of the samples are known, the 
rest of the data will not provide any information about the target because 
any sensor in between must be a detecting sensor. For the sake of notational 
simplicity, let us denote zx^^ and z-ou hy x and y respectively in this sub¬ 
section. If we consider any estimator of z^^ as Fp = the Rao-Blackwell 

theorem indicates that h{x,y) = E[v{Zd) \ {x,y)] would have a smaller or 
equal variance to that of ^(Zd) where E represents the expectation of the 
random variable HU. Thus, we can restrict our search for an MVU estimator 
to the functions of {x,y). 

Although (x, y) are a minimal sufficient statistics for this problem, a com¬ 
plete sufficient statistics does not exist iSlIlZ!. For the case of a uniform 
distribution, |12] proved that no non-constant parametric function of Zt (in¬ 
cluding Zt) can admit a Uniform Minimum Variance Unbiased estimator 
(UMVUE). |l8] provided another proof for the nonexistence of UMVUE for 
this problem which is true even if n = 1 (estimators built based on one obser¬ 
vation). However, in location estimation the estimators have characteristics 
that limit the function space and a Minimum Variance Unbiased estimator 
may exist among those functions. For example, in location estimation (and 
perhaps many other estimators of physical phenomenon) it is reasonable to 
assume that estimators are invariant under Euclidean motion i.e. we expect 
that if the observations shift, then the estimate of the phenomenon shifts 
along with the observations as illustrated in Figure |C.l . Therefore, if we 
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estimate zt by h(x, y), we should have 

h{x + d,y + d) = h{x, y) + d. (1) 

Let us dehne the possible target region for as T(Zd) = [y — R,x + R] 
as illustrated by hashed interval in Figur ^C^ From the dehnition, it is clear 
that T(Zd) = T{x,y). An observation {x,y) implies that zt G T{x,y). On 
the other hand, any observation in which z^ G [y — R, x + R] is a possible 
observation because this along with x < y guarantee that \x — zt:\ < R 
and \y — zt:\ < R. It can be shown that the T{x,y) interval and {x,y) are 
related by a one to one mapping. Thus, in the search for an MVU function 
it is possible to search on functions of T {x, y) instead. In other words, the 
estimator can be limited to. 


giT{x,y)). 

Therefore, the invariance under Euclidean motion suggests that 
g{T{x + d,y + d)) = g{T{x,y)) + d. 


( 2 ) 

(3) 


Now consider a sub-partition of random space that has n detecting sensor. 
The mean squared error for this sub-partition becomes. 


E[\zt 



r*2T+^ /y _ ^ 


' z^—R J Z'j^—R 


n—2 



\ 9 {'T{x,y)) - Zt 


lx<ydydx. (4) 


where lx is the indicator function of X, the binomial coefficient compensates 
for the fact that any two of the n sensors can be selected as the maximum and 
minimum values of the data; and is the probability that all other 

sensors are located inside [x^y]. Now, with a change of variable t = y — x, 
removing y and considering that t G [0,2i?], we have, 

rz-r+R r'^R / t / 1 \‘^ 9 

E[\zt-Zt\ ] = n{n-l) (^-j \g{r{x,x + t)) - zt\ 

'^ZT^[t^x—R,R^x]dtdx. (5) 

The assumption regarding invariance of g to Euclidean motion results in 
g{T{0+x,t+x)) = g{TiO,t))+x. In addition, l^^(.[t+x-R,R+x] = lzT-xe[t-R,R]- 
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Thus, we have, 


/■zt^+R |•2R / f \n 2 / ]^ \ 2 

E[\zt-Zt\ ] =n{n-l) (^-J \g{r{0,t)) + x - Zt\ 

'^ZT—xG[t—R,R]dtdx. ( 6 ) 

With another change of variable t' = — x and eliminating x and we 

have, 

E[\zt- ZT\]=n{n-l) (-j |^(r(0, t)) - t'| 

dtdt'. (7) 


Moreover, because the interval of the two integrals are independent, we can 
change the order of the integrals and obtain. 


f2R / f \ n—2 / 1 \ ^ p+R „ 

EWzt- z^\]=n{n-l) (-j l^-J jg(r(0,t))-t'l lt'e[t-R,R] 

dt'dt. ( 8 ) 

Furthermore, because f > 0 we know that [t — R, R] C [—R, R], thus. 


p2R / f \ n—2 / 1 \ ^ pK 

E[\zt-zt\] = n{n-l)J^ (-j f-j J^^Jg{r{0,t))-t'\dt' 


rR 

It-R 


dt. 

(9) 


Substituting T(0, t) by its value [t — R, R] we have. 


E[\zt — zt^] 



r \gi\t-R, R])-t'\^dt' 

Jt-R 


dt. 

( 10 ) 


On the other hand, g{[t — R, i?]) does not depend on t' and we know that if 
it is selected as the center of gravity of the inner integration interval i.e. 

g(\t - R, fi|) = CG{\t - R,R]) = (n) 

the inner integral will always be minimized [ISHSI] • Hence because the outer 
integral integrand is always positive, the entire expectation will be minimized. 
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Furthermore, similar to derivation presented in (10), we can show that 


E[zt: — z-r] 



r {g{[t-R, R])-t')dt^ 
Jt-R 


dt. 

( 12 ) 


Thus, it can be verihed that g{[t — R, i?]) = CG{[t — R, i?]) = * makes 
(12) zero i.e. (11) is unbiased. Considering that the following are equivalent 


g{[t + X — R, X + R\) = g{[t — R, R]) + x = CG{[x + t — R, x + R]) (13) 


Then, g(T{x,y)) = GG{T{x,y)) will be the minimum variance unbiased 
estimator among all estimators that are invariant to Euclidean motion. 


3.2. Unknown Detection Radius 

In this sub-section, we will consider the problem when R is unknown 
which means that [z^, R] are both parameters of estimation now. In this 
case, {zbltZbu) iiof is rninimal sufficient statistics but also is complete 
[inillZ]. Now consider a mean estimator u(Zd) = ^ According to 


to Rao-Blackwell theorem 


Zt — E [v{Z-[))\{z-DL, ^Bu)] — 


^DL+^Df/ 


would 


be the MVU estimator which is the same result as has been derived in Section 


3.1 but without any constraint on function space of the estimator and for a 


different reason. 


4. On the Existence of an MVU Estimator in the Non-Censored 
Scheme 

In this scenario both detecting and non-detecting sensors report their 
locations and observations to FC. Therefore, the problem can be viewed as 
estimation of the center of a uniform distribution over [zt—R, zt+R] given its 
samples and samples of another uniform distribution over the complemented 
interval [Ai, A 2 ] — [zt — R, zt + R]. We dehne the following, 

Zbl = min ZBi 

*SiSd 

zbu = m&xzBi 

iScSd 

^ max z^bj 

JSOnD & Zjq£,j<ZDL 

^NDu = . ^ min zndu (14) 

& ^NDj>2:D(7 










With the assumption that at least one sensor is a detecting sensor, it’s 
straight forward to show that 2 ;ndl summarizes the information of all non¬ 
detecting sensors located at ^ < zdl, and summarizes the information 
of all non-detecting sensors located aX z > z-ojj. In addition, the combination 
of zdi and zjm summarizes the information of all detecting sensors. 


4 . 1 . Known Detection Radius 

In this scenario, the possible target region can be dehned as 

7~(Zd, Znd) = [zbu ~ R-, ^DL + i?] — ([^^NDL — R-, + -R] U 

[2^ND(7 ~ R: + -R]) (15) 

The following conditions categorize the possible situation of zdl, zdu, z-^bl 
and 2 ;nd ;7 with respect to each other as illustrated in Figures C.3[a)-(d). 


*) ~ R) ^ (^NDL + .R) & i^BL + -R) ^ (^NDf/ ~ R) 

ii) {zbu — R) < (^^NDL + i?) & {zbl + R) > ~ R) 

Hi) {zbu — R) < (2:ndl + i?) & {zbl + -R) < — R) 

iv) {zbu - R)> (^NDL + R) ^{zblR R) > - R) 

Depending on which of these conditions are true, T(Zd, Znd) could be calcu¬ 
lated based on only two of the elements of {zbd zbui z^bli ^ndc/} as follows 


7~(Zd, Znd) 


[zbu ~ R-, Zb lR R\ 
[^^NDL + R-, ^NDi7 — .R] 

[^NDL + .R, ^DL + .R] 

[zbU ~ Ri ^NDC/ ~ -R] 


if condition{i) 
if condition{ii) 
if condition{iii) 
if condition{iv) 


(17) 


Let us represent the beginning and end of the interval T(Zd,Znd) by 
Ti and T 2 . 7~(Zd, Znd) = [DjD] incorporates all information regarding 

the target location. Moreover, an observation is valid if and only if Z't G 
T(Zd, Znd)- It can be shown that, T(Zd, Znd), is a minimal sufficient statis¬ 
tics for estimation of zt although it is not complete. 


^To show the incompleteness we can equivalently consider the minimal sufficient statis- 
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Moreover, the probability density function that certain detecting sensors 
are located at Zd and non-detecting sensors are located at Znd, can be 
dehned as 

/ 1 \ 

/(Zd, Znd) = IzteTCZd.Znd) (18) 

Because the four conditions in 0 can not be simultaneously true for any 
observation, we can partition the observation space based upon which condi¬ 
tion is met. We will proceed and show that g = CG (T(Zd, Znd)) is unbiased 
and minimizes the MSE on each of these partitions under the condition that 
the deployment region is sufficiently large and the estimators are invariant 
to Euclidean motion. Therefore, g = CG (T(Zd, Znd)) is an MVU estimator 
over the entire random space. 

We will provide the proof for one of the partitions below. The details 
of the other three are similar and have been provided in appendix. We will 
partition each sub-space of random space further by considering the number 
of detecting sensors and non-detecting sensors. 

4 .I.I. proof for condition Hi 

This condition means that T(Zd, Znd) = [-^ndl + R) ^dl + R] ■ In this 
sub-section, let us denote ^ndl and zdl by x and y respectively. Hence, 
the possible target region would be T(Zd, Znd) = 7~(a;, y) = [x + R,y + R], 
Moreover, we have that 


Zdu < ^NDL + 2i? 

^ND{/ > ^DL + 2i? 

zdu ^ (19) 

Therefore, all detecting sensors should be located in interval [y, x + 2R] and 
all non-detecting sensors should be located in interval [Hi, x] or [y -|- 2R, A 2 ]. 
An {x,y) observation in this case implies that zt G T{x,y). Note that in 
this case y < z^u < x + 2R along with zt G [x + R, y + R] guarantees that 
Zt — 3R < X < Zt — R and [zt — R) < y < [zt + R)- 


tics as (ri, (t 2 — ri)). In this way, (r 2 — ti) is an ancillary statistics whose PDF does not 
depend on zt but in combination with ti provides information about zt- Therefore, the 
minimal sufficient statistics fail to rid itself from ancillary part of data which might be 
considered as an indication that it is not complete [46) . 
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As explained in Section 3.1 


we can restrict onr search for an MVU es¬ 
timator to fnnctions of the minimal snfficient statistics i.e. g{'T{{x,y))) = 
q{\x + R,y + R\). The invariance to Euclidean motions result in qiTix -\- 
d,y + d)) = g{T{x,y)) + d. 

Now consider a sub-partition of partition iii of random space that has p 
non-detecting sensors and n detecting sensors. The mean square error for 
this sub-partition would be, 


E[\z\ 


'L-2R 


n' 


{y 


1 

L 

X] 


2 nZT-R rZT+R / 2R — (y — x) 

RZ'ySR \ L 

P-1 


n—1 


\ 9 {T{x,y)) - zt\ ly-x< 2 Rdydx ( 20 ) 


where jg probability that all other detecting sensors are 

located inside [y^x -|- 2R\ and jg ^pg probability of all other 

non-detecting sensors will be located in [Ai,a:] U [?/ -|- 2i?, A 2 ] interval. Now, 
with a change of variable t = y — x, removing y and considering that t G 
[0, 2R\. We have. 



\g{T{x,x + t)) - zxf l^^(zi^+R^^+t+R]dtdx (21) 

where lzr^^[x+R,x+t+R] and t G [0, 2R\ guarantees that y will be located in [z^ — 
R^zx + R]- The assumption regarding invariance of g to Euclidean motion 
result in g{T{0 + x,t + x)) = g{T{0,t)) x. In addition, = 

l^T-a;e[R,t+R]- Thus, we have. 


E[\zt — zxf] 


np rz^-R r 2 R /2R - ty ^ /L - 2R - ty-^ 

Jzt:- 3 R Jo \ L J \ L ) 
\giT{ 0 , t))+x- zxf lz^-xe[R,t+R]dtdx 


( 22 ) 


With another change of variable t' = zx — x and eliminating x and zx we 
have. 


E\\zx — z^] 


np r5R r2R /2R - ty-^ /L-2R- ty-^ 

iR Jr Jo V L ) V L ) 

\g{T{ 0 R)) - lt'(z[R^t+R]dtdt' 


( 23 ) 
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Moreover, because the interval of the two integrals are independent, we can 
change the order of the integrals and obtain. 


E\\z-£ — z^f] 


np r2R rSR /2R - /L-2R- ty-^ 

1? Jo Jr V L J V L J 

\g{T{0,t)) - t'\‘^ lt'elR,t+R]dt'dt 
np /2R-ty-WL-2R-ty-^ r^R 

1? Jo V L ) V L ; Jr 

\g{T{0,t)) - t'f lt'eiR,t+R]dt'dt 


(24) 


Furthermore, because t G [0, 2R] we know that [i?, t + R] G [-R, 3R], thus. 


j?\\ - ih npr^Rf2R-t 

El\zr - zr\ ] ^ I 


n—1 


L-2R-ty-^ fi+R 


-) 


IR 


\g(T{0, t)) — dt'dt 

Substituting T(0, t) by its value [R, R + t] we have, 

^L-2R-t^P-^ 


77 n - ih npr^Rf2R-t 

El\zr - z-rl ] = 


L 

■t+R 


rt+n 

/ \g{[R, R + t])-t'fdt' 

^ R 


dt 


( 26 ) 


( 26 ) 


On the other hand, g{[R, R + t]) does not depend on t' and we know that if 
it is selected as the center of gravity of the inner integration interval i.e. 


g{\R, t + R])= CG{\R, t + R]) = 


t + R + R 


(27) 


the inner integral will always be minimized |19H5T] . Hence because the outer 
integral integrand is always positive, the entire expectation will be minimized. 
Furthermore, similar to derivation presented in (26), we can show that 

np GR(2R-G^-^ 


E[+t — zt] — 


L2 


L 

-t+R 


L-2R-ty-^ 




/ {g{[R, R + t])-t')dt' 

JR 


dt 


( 28 ) 


Therefore, it can also be verihed that g{[R., R + t]) = CG{[R, i? + f]) = 


^ makes (28) 


zero i.e. 


the estimator is unbiased. Considering that the 


12 



















followings are equivalent 

g{[x + R, X + R + t]) = g{[R, R + t]) + x = CG{[R, R + t]) + x 

= CG{[x + R, x + R + t]) (29) 

Then, g{T{x,y)) = GG{T{x,y)) will be the minimum variance unbiased 
estimator among all estimators that are invariant to Euclidean motion. 


4-2. Unknown Detection Radius 

In this scenario, we will reconsider the results when R is itself a target 
of estimation. Similar to Section 13.21 the most informative sensor locations 
are 2 ;dl, z-du^ ^ndl and However since R is not known, dehnition of 

possible target region is not as straight forward as before. Infact in this case, 
R can admit a range of values if it meets following conditions. 


R > 


Zbu + Zt>l 


R < 


2^NDL + 


(30) 

(31) 


Meeting the first condition means that interval [zbu ~ R-, zbl + -R] in (15) is 
not empty and the second condition guarantees that the balls around non¬ 
detecting sensors zndl and ^ndc/ will not exclude all the interval from being 
possible target region. 

Now, consider a valid observation where all zbl, zbu, 2^ndl and 2 ;nd (7 
exist. Depends on what is the actual value of parameter R, a different condi¬ 
tions in (16) is met and T and g functions would be different. Let us dehne 
gn (Zd, Znd) as follows: 


dR (Zd, Znd) — 


2 

^NDL+l^NDf/ 

2 

j 

^DrZ+l^NDr/ _j_ 


if 

if 

if 

if 


ZDr/-^NDL 


k{R< 


{R > 
{R > 

{R < 


z-pu—z-tm L 
2 

zpu—z^u £ 
2 

^Dr/~^NDL 

2 


^{R> 
& (i? < 


ZlSIPU—ZPL 

2 

Z^PlJ — Zp 
2 


(32) 


& (fl > amipDi) 


Let us consider two different values for R, Ri and R 2 such that each of them 


enable a different condition in (32). According to Section 4.1 gR^ (Zd,Znd) 
and gR.^ (Zd, Znd) would be the MVU estimator when R = Ri and R = R 2 


respectively but they might have different values as illustrated in Figure C.4 


Thus, an MVU estimator does not exist for R unknown case even if we restrict 
the function space to those Invariant to Euclidean motion. 
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5. Conclusion 

We considered the problem of localization with binary ideal RSS measnre- 
ments in this paper. Finding a Minimnm Variance Unbiased estimator for 
this scenario is important because it establishes a lower bound for all other 
estimators in the presence of uncertainty. We showed that even though in the 
theoretical statistics literature a version of the problem has been categorized 
as a problem that does not admit any MVU estimator, the result might be 
different when other characteristics of the estimators are considered. We add 
the assumption of non-invariance of the estimator under Euclidean motion 
and show that when the detection radius is known the MVU estimator does 
exist for both censored and non-censored scenarios under this constraint, 
whereas it still does not exist for the non-censored scenario when the detec¬ 
tion radius is unknown. We believe that this constraint is critical for many 
other applied estimators, thus care is necessary when applying pure theoret¬ 
ical statistical results to estimation of physical phenomenon. Hence, it may 
justify the necessity of re-evaluation of some of these results this assumption 
(or perhaps other ones too) regarding the functional space. 

Appendix A. proof for condition i 

This condition means that T(Zd,Znd) = —-R, +-R]- In this 

sub-section, let us denote and z-qu by x and y respectively. The possible 
target region would be T(Zd,Znd) = T{x,y) = [y — R,x + R\. Moreover, 
under this condition we have that 


^NDL < zj:)u — 2R 

2^ND(7 ^ ^dl + (A.l) 


Therefore, all non-detecting sensors should be located in interval [Ai, z-du ~ 
2R] or [zdl + 2i?, A 2 ]. An observation {x,y) is valid if and only if z^ € 
[y — R,x + R]. 

As explained in Section 3T, we can restrict our search for an MVU es¬ 
timator to functions of the minimal sufficient statistics i.e. g{T{{x,y))) = 
g{[y — R,x + R]). The invariance to Euclidean motion assumption suggests 
that g{Tix + d,y + d))= g{T{x,y)) + d. 

Now consider a sub-partition of random space that has p non-detecting 
sensors and n detecting sensors. Let us assume that n > 1. We will examine 
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the result when n = 1 at the end of this sub-section separately. We can write 
the Mean Squared Error as following: 


E[\z^ — Z'rf' 


2 (n\ fy - f L - AR +y - x\^ 

^V2; Jzt-R Jz-r-R V L ) V L ) 

\giT{x,y)) - zt\^ lx<ydyd^A.2) 


where the binomial coefficient compensates for the fact that any two of the 
n detecting sensors can be selected as the maximum and minimum values of 
the locations of the detecting sensors and is the probability that 

the other detecting sensors are located inside [x,?/]; and ig 

probability of all non-detecting sensors will be located in [Ai, y — 2R\ U [x -f 
2R, ^ 2 ] interval. Now, with a change of variable t = y — x, removing y and 
considering that t G [0,2i?]. We have, 


E\\z^ — ZT:f‘ 


n(n — 1) r^T+R ^ f L — 4:R + t V 

lA Jzr-R Jo ^l) V L ) 

\g{T{x,x + t)) - zt^ lzT,&[t+x-R,R+x]dtdx 


(A.3) 


where lzr^^\t+x-R,R+x\ guarantees that x and x + t will be located in the de¬ 
tecting area. The assumption regarding invariance of g to Euclidean motion 
result in g{T{^ + xR + x)) = g{T{^R)) + x. In addition, l^^<z[t+x-R,R+x] = 
Thus, we have. 


E\\z^ — z^'^ 


n{n — 1) r‘2R / t W ^ /L — AR + t\^ 

Z2 Jzr-R Jo ^LJ V L ) 

\g{T{0, t))+x- zt^ l^^_^(zit-R^R]dtdx 


(A.4) 


With another change of variable t' = z^ — x and eliminating x and z^ we 
have. 


E[\zt 


i 2 t n(n 
= 


(n — 1) f+^ ( t 


L2 


l-R Jo 


n—2 


L-AR + ty 


lt'&[t-R,R]dtdi^A.5) 
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Moreover, because the interval of the two integrals are independent, we can 
change the order of the integrals and obtain. 


E[\z^ — zt^] 


n{n-l) /L-AR + ty r+R 

Jo yi) V L J J-R 

\g{T(l),t)) - t'f 


(A.6) 


Furthermore, because 0 < t < 2R we know that [t — R, R] C [—R, R], thus, 

n{n-l) fL-AR + t^P 


Euzt — zt| — 


L2 


Jt-R 


dt 


(A.7) 


Substituting T(0, t) by its value [t — R, R] we have. 




^-2 ^L-AR + ty 


r \g{\t-R, R])-t'\^dt' 
Jt-R 


dt 


(A.8) 


On the other hand, g{[t — R, i?]) does not depend on t' and we know that if 
it is selected as the center of gravity of the inner integration interval i.e. 


g{[t-R, R]) = CGi[t-R,R]) = 


(A.9) 

the inner integral will always be minimized [ISHSI] • Hence because the outer 
integral integrand is always positive, the entire expectation will be minimized. 
Furthermore, similar to derivation presented in (A.8), we can show that 

nin — 1) ft'' 


E[z^ — zt\ — 


L2 


rR 


' t—R 


{g{[t-R, R]) 


L-AR + t 
' L 

f) dt' 


dt 


(A. 10) 


Therefore, it also ca n be verihed that g{[t — R, i?]) = CG{[t — R,R\) = 
f-fi+R (A.10) zero i.e. the estimator is unbiased. Considering that 

the followings are equivalent 


g{[t + X — R, X + i?]) = g{\t — R, R]) + x = GG{[t — R, i?]) + x 

= GG{[x + t-R,x + R]) (A.ll) 
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Then, g{T{x,y)) = CG(T{x,y)) will be the minimum variance unbiased 
estimator among all estimators that are invariant to Euclidean motion. If 
n = 1, the Mean Squared Error would become 

„ 1 r^T+R /L - 4:R\P 

E[\zt - zt\ ] = - J ^ — - —j \g{T{x,x)) - zt\ dx (A.12) 

with a change of variable t' = z^—x and similar effort as previously presented, 
we will get that, 

o 1 r+^/L-4:R\^ 9 

E[\z^ - ^ j \9im 0)) - tf dt' (A.13) 

Therefore, 5'(T(0, 0)) = 0 make the MSE minimized and therefore 5 '(T(a;, a;)) = 
X will be the MVU estimator. 

Appendix B. proof for condition ii 

This condition means that T(Zd, Z^d) = [^^ndl + R, — .R]- In this 
sub-section, let us denote ^^ndl and ^^NOf/ by x and y respectively. The 
possible target region would be T(Zd,Znd) = 'T{i.x,y)) = [x + R,y-R]. 
Moreover, 


Zdu < ^NDL + 2i? 

Z-DL> ^NDC/ — 2i? (B.l) 


Therefore, all detecting sensors should be located in interval {y — 2R,x + 2R) 
and all non-detecting sensors should be located in interval [Ai,a;] or [y, A 2 ]. 
In this case, an observation {x,y) implies that z-r G [x + R,y — R\. Note 
that in order for [y — 2R, x + 2R\ to be a valid interval it is implicit that 
0 < y — X < AR. 

we can restrict our search for an MVU es- 


As explained in Section 3.1 


timator to functions of the minimal sufficient statistics i.e. g{T{{x,y))) = 
g{[x + R,y — R]). The invariance to Euclidean motion assumption suggests 
that g{T{x + d,y + d)) = g{T{x, y)) + d. 

Now consider a sub-partition of random space that has p non-detecting 
sensors and n detecting sensors. We can write the Mean Square Error as 
following: 
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-E'[|iT — zt \'^ 


2 /p\ rzT-R .zt+sr /L-{y-x)Y~^ 

£^[2) Jz^-sRjzr+R V L ; V L ; 

\g(T{x,y)) - zY ly-x<4RdydiB.2) 


where the binomial coefficient compensates for the fact that any two of the 
p non-detecting sensors can be 2 :ndl and ;2 ndc/- Moreover, is 

the probability that all detecting sensors are located inside [y — 2R, x + 2R] 
and is the probability of all other non-detecting sensors will be 

located in [Ai, x] U [y, A 2 ] interval. Now, with a change of variable t = y — x, 
removing y and considering that t G [2i?, 4i?], we have. 


E \z'j' — Z'j' 


p{p-i) y-R yR - ty y - ty-^ 

Jzt-3R J 2 R V L ) \ L ) 
\g{T{x,x + t)) - zY l^^(.[^+R^^+t_R]dtdx. 


(B.3) 


where lzr^^\x+R,x+t-R\ guarantees that x and x + t will be located in the non¬ 
detecting area. The assumption regarding invariance of g to Euclidean mo¬ 
tion result in 5f(T(0-|-a;, t-l-a;)) = g{T{QR))+x. In addition, lzrj,^[x+R,x+t-R] = 
l^^_x(i[R,t-R\- Thus, we have. 


BN* _ ,^ 1^1 = p(p - 1 ) /■"-* r (iR -1 


L2 


' — ^ 2R 


L-t 


p-2 


\g{T{0,t)) + X - zY l,^_xGiR,t-R]dtdx. (B.4) 

With another change of variable t' = z^ — x and eliminating x and z^ we 
have, 

Bii* - .,1=1 = p(p -1) - n Vi - 


Jr J 2 R \ L J \ L 

\g{T{0,t)) - t'l"^ lt'(z[R^t-R]dtdt' (B.5) 

Moreover, because the interval of the two integrals are independent, we can 
change the order of the integrals and obtain. 


J 2 R 


IR 


\g(T{0,t)) -t'f lt'eiR,t-R]dt'dt (B.6) 
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Furthermore, because t G [2R,4:R] we know that [R, t — R] C. [-R, 3-R], 
thus, 


pn ^ |2n pip - 1) f4:R-t 

E[\zt - = 


L2 


'2R 
rt—R 

Jr 


\9iTm) 


L 

t'\^ dt' 


dt 


(B.7) 


Substituting T(0,t) by its value [R,t — R] we have, 

pip-1) /AR-t\^ fL 


E Ft — — 


L2 


l2R 


P-2 


f ^\gi[R, t-R])-t'\^ dt' 

J R 


dt 


(B.8) 


On the other hand, giiR, t — -R]) does not depend on t' and we know that if 
it is selected as the center of gravity of the inner integration interval i.e. 

gi[R, t-R]) = CGi[R, t - R]) = (B.9) 

the inner integral will always be minimized [inHSI]. Hence, because the 
outer integral integrand is always positiv^ the entire expectation will be 
minimized. 

Furthermore, similar to derivation presented in (B.8), we can show that 
E\z^ — '^tJ — 


L2 


I2R 
rt—R 


f {g{[R, t-R])-t')dt' 

Jr 


dt 


(B.IO) 


_ t-R+R 


makes 


It also can be verihed that g{[R, t — i?]) = CG{[R, t — i?]) = 

(B.IO) zero i.e. the estimator is unbiased. Considering that the followings 
are equivalent 


gi[x + R, X + t — R]) = gi[R, t — i?]) + x = GGi[R, t — R\) + x 

= GGi[x + R, x + t-R]) (B.ll) 

Then, giTix,y)) = GGiTix,y)) will be the minimum variance unbiased 
estimator among all estimators that are invariant to Euclidean motion. 


^Note that L — t > 0 because of the assumptions that the deployment is so large that 
[zt — 3i?, zt + 3i?] C [Ai, A 2 ] 
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Appendix C. proof for condition iv 

Conditions iv) would be similar to conditions iii) because of symmetry. 

Therefore, we will not repeat the proof here for the sake of space. 
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Figure C.l: Non-invariance to Euclidian motion: zt = zt + d 
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Figure C.2: Possible target region ilustration for censored scenario, ri = 
Zdu - R and T 2 = z^l + R. 
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Figure C.3: The four possible relative locations of ^^DLXDc/XNDLXNDn in 
non-censored scenario and the effect in ti and T 2 determination. 
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27 



